Vortex states of the E u model for Sr 2 Ru0 4 
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Based on the Ginzburg-Landau functional of E u symmetry presented by Agterberg, vortex states 
of Sr2Ru04 are studied in detail over H a i J$ H < H C 2 by using the Landau-level expansion method. 
For the field in the basal plane, it is found that (i) the second superconducting transition should 
be present irrespective of the field direction; (ii) below this transition, a characteristic double-peak 
structure may develop in the magnetic-field distribution; (iii) a third transition may occur between 
two different vortex states. It is also found that, when the field is along the c axis, the square vortex 
lattice may deform through a second-order transition into a rectangular one as the field is lowered 
from H C 2. These predictions will be helpful in establishing the E u model for S1-2RUO4. 



Active studies have been performed on superconduct- 
ing Sr2Ru04 Q where another unconventional pairing 
may be realized [Q. A possible candidate for its sym- 
metry is the E u model with two-fold degeneracy 0, as 
indicated by various experiments [^~[To|| . However, fur- 
ther experiments seem being required before establishing 
its validity for Sr 2 Ru04. In this respect, the vortex states 
may provide clear and indisputable tests for the p-wave 
hypothesis. The present paper provides a detailed theo- 
retical description of them which will be helpful towards 
that purpose. Clarifying the basic features of the two- 
component model, which has not been performed com- 
pletely, will also be useful for the experiments of UPt3. 

The vortex states of the E u model for S^RuCU 
have been studied theoretically in a series of papers by 
Agterberg et al. |ll| ^3|. Based on the two-component 
Ginzburg-Landau (GL) functional and following essen- 
tially Abrikosov's method (li) which is effective near the 
upper (H C 2) and lower (H c %) critical fields, they have 
provided several important predictions. Especially note- 
worthy among them are: existence of the second transi- 
tion for H _L c similar to that observed in UPt3 |l5],|l6| ! 
several orbital-dependent phenomena helpful in identify- 
ing which band is mainly relevant; stabilization of the 
square vortex lattice for H | c. An observation of the 
square lattice has been reported by Riseman et al. [ pT| . 

With these results, this paper focuses on the following: 
(i) The properties of the intermediate fields, in particular 
those below the second transition for H_Lc, remain to be 
clarified. We will treat the whole range Hci H H C 2 in 
a unified way, describe possible changes of experimentally 
detectable properties as a function of the field strength, 
and draw characteristic features in low fields, (ii) Con- 
sidered have been the cases where the field is along the 
high-symmetry axes. It is still not clear whether or not 
the second transition for H_Lc persists for arbitrary field 
directions in the ab plane, because the term |?7i| 2 |?72| 2 in 
the GL functional [see Eq. (Q) below] generally causes the 
first- and third-order mixing. We will study those gen- 
eral cases to establish the existence of the second tran- 



sition, (iii) Agterberg introduced several assumptions in 
the parameters used to minimize the free energy. We will 
perform the minimization without such assumptions. 

The goals (i)-(iii) may seem rather formidable, but 
they can be achieved with the Landau-level expansion 
method jL8|. When applied to the s-wave pairing, it 
successfully reproduced the properties of the whole re- 
gion H c i < H < H C 2 quite efficiently for an arbitrary k. 
Compared with the direct minimization procedure in real 
space [|9| , the method has a couple of advantages that (i) 
it is far more efficient and (ii) one can enumerate possi- 
ble second-order transitions rather easily, hence enabling 
us to establish the phase diagram of various multi-order- 
parameter systems. This is the first time where it is 
applied to an multi-order-parameter system so that this 
paper also has some methodological importance. 

The GL free-energy density adopted by Agterberg is 
given by O 
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where the same notations are used here. This simplified 
free energy has the advantage that there are only two 
parameters in it whose values can be extracted from ex- 



1-37 
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the latter 



periments, i.e. Ki = H C 2 /v2H c and 
being related to the H C 2 anisotropy in the ab plane as 
H c2 (a)/H c2 (a + b) = 1=£ 0. The value « x = 31 (1.2) 
for H_Lc (H || c) will be used throughout pQ], whereas 
v is left as a parameter. A recent observation of the H C 2 
anisotropy suggests that v is positive and - 0.01 10. 

We sketch the method to find the minimum for an arbi- 
trary field strength |0| . Let us fix the mean flux density 
B = (B sin 6* cos ip, B sin sintp, B cos 6*) rather than the 
external field H, and express h = B + h where the spatial 
average of h vanishes by definition. We then transform 
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where <^>, 0' and L are conveniently chosen as <f> — 
tan" 1 [27 tan 2^/(1 - 7)], 0' = L 2 cos9, and L = 

{(l+7-/)/[(l+7-/)cos 2 + 2 K5 sin 2 0]} 1/4 with / = 

1/2 

[(1 — 7) 2 cos 2 2iy9 + 47 2 sin 2 293] . Assuming uniformity 
along z' direction, we then expand rj'(r') and h (r') as 



V '(r>) 
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h(r') =z'J2 h * exp(iK • r') , 



(4) 
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where is the system volume, ipNq denotes an eigenstate 
of the magnetic translation group in the flux density B 
with the Landau-level index N and the magnetic Bloch 
vector q, and K is the reciprocal lattice vector of the 
vortex lattice. The explicit expression of ^jvq( r ') for the 
spacial case where one of the unit vectors of the vortex 
lattice, a 2 , lies along the y' axis is given by 

A/f/2 

'0A'q( r ')= c^ q y' (y ' +0 - 5q ^^ +nai ^^ v ' + ' !: ''~ - 5nai «'^ /l ^ 
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with A/" 2 the number of flux quanta in the system, l c 
denoting -h= of the magnetic length, and a\ x * (ai y >) the 

x' [y') component of another unit vector ai Eq]. We 
also consider the counterclockwise rotation of ai and a 2 
around the z' axis by the angle ip' L . Substituting Eqs. 
@"@ m t° Eq. (0) and integrating over the volume, we 
obtain the free energy per unit volume as 



F[{c N J , {h K }, B, p, V, <p' L ] = i //[■//( r' ). 1, ( , ' ). B\,l 
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where p = |ai|/|a 2 | and d = cos 



-1 ai-a 2 



(6) 

This F is a 



. |ai11 ? 21 ' 

desired functional which can be minimized rather easily 
using one of the standard minimization algorithms |22| . 
Due to the periodicity of the vortex lattice, we only have 
to perform the integration over a unit cell. The external 
field H is then determined through the thermodynamic 
relation (H = in the present units). In numerical 
calculations we have cut the series in Eqs. (0) and (H) 



at some N c and |K C |, respectively, thereby obtaining a 
variational estimate of the free energy. The convergence 
can be checked by increasing N c and |K C |. The choice 
N c ~ 12 and K c ~(the third smallest) has been checked 
to provide correct identification of the free-energy mini- 
mum with the relative accuracy of 10 -6 for B / £f C 2 >0.1. 
Though not presented here, preliminary calculations re- 
veal that the method is also effective for #^0, 

The functional F has another advantage that one may 
enumerate possible transitions in the vortex states of 
multi-order-parameter systems. Much attention has been 
focused on this subject in connection with the observed 
phase diagram of superconducting UPt3 |23|]. No com- 
plete analysis has appeared yet, however, and the use of 
F will be quite helpful for that purpose. The features of 
the s-wave lattice can be summarized as follows [Q: (a) 
a sing le q in Eq. (|) suffices to describe it with a choice 
of q corresponding to the broken translational symmetry 
of the lattice; (b) the hexagonal (square) lattice is made 
up of N = 6n (4n) Landau levels (n: integer) |24|,|ll|; (c) 
more general structures can be described with N = 2n 
levels, odd N's never mixing up since those bases have 
finite amplitude at the core sites; (d) the expansion co- 
efficients cjvq can be chosen real for the hexagonal and 
square lattices. With these results on the conventional 
lattice, the following second-order transitions are possi- 
ble in multi-component systems: (i) deformation of the 
hexagonal or square lattice which accompanies entry of 
new iV's as well as complex numbers in the expansion 
coefficients; (ii) mixing of another wave number q 2 satis- 
fying q 2 -qi =K/2 Q2|; (iii) entry of odd N's. Though 
not complete, this consideration will be sufficient below. 




FIG. 1. Transition lines for H || a as a function of the 
anisotropy parameter v. The closed (open) circles denote the 
zeros of r\\ The inset plots the angle <p' L at the I< — >II 

transition as a function of v. 
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FIG. 2. Spatial variation of l z = 
state II for v = 0.077 and B/H c2 
l z ~ ±1 correspond to the "bulk" 
corresponding |»j(r)|. 
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FIG. 4. The magnetic-field distribution P(h) in the state 
II for B/H c2 = 0.45, 0.35, 0.25, and 0.15. A characteristic 
double-peak structure develops as the field is decreased. 
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FIG. 3. A comparison of |r/(r)| between the states II and 
III with v = 0.077 and B/H c2 = 0.45, 0.25. The amplitude 
| ?7 (r) | is finite everywhere in the state II, which is brought 
about at the expense of the variation in l z ; see Fig. [| 

We now present the results for H _L c. Figure 1 shows 
the transition lines for H 1 1 a (6 = ^ ; (p = 0) as a function of 
the anisotropy parameter v\ the one given as a function 
of 7 = has qualitatively the same structure, with 
7 = and 1 respectively corresponding to v = 1 and — 1 . 
As already pointed out by Agterberg [jll| , there are three 
possible vortex states: the high-field region I where a hex- 



agonal lattice is stable with 772 = 0; the region II where 
rj2 becomes finite with c\2 — Qi equal to half the unit vec- 
tor bi of the reciprocal lattice, i.e. the vortex lattice is 
coreless with |tj| finite everywhere; the region III where a 
deformed conventional lattice with 77 2 7^ is stable. In ad- 
dition, Fig. 1 includes the following new results: (i) a full 
minimization with respect to <p' L clarifies that the I<->II 
transition is continuous as a function of v (see the inset); 
(ii) high-precision calculations in the low-field region re- 
veal that, as the field is lowered, the coreless state II is 
replaced via a first-order transition by the state III with 
cores. The reason for (ii) can be realized by looking at the 
variation of l z = (r/xr)* )-z/2i\r]i\\r}2 \ which is proportional 
to the magnitude of the orbital angular momentum along 
z. As seen in Fig. 2 calculated for v = 0.077 (7 = 0.3) and 
B / H c2 = 0.25, one of the bulk states l z = H is alternately 
realized in II, and there necessarily exist lines of "defects" 
where l z vanishes. Compared with III where \rj\ vanishes 
at points, the state II is thus energetically unfavorable at 
low fields. It can however be stabilized at intermediate 
fields by making \rj\ more uniform. Figure 3 plots \t](t)\ 
for v — 0.077, showing how the differences between II 
and III develop as B/H C 2 is decreased. In fact, only a 
deformation of the lattice occurs in III, whereas a layered 
structure also shows up in II with |?7(r)| becoming more 
and more uniform. This rather drastic change in II can 
be detected by measuring the magnetic-field distribution 
P(h) = yj S[h — h x (r)} d 3 r. As seen in Fig. 4, the sin- 
gle peak at B/H C 2 = 045 splits and one of them moves 
towards the high-field end, which originates from the de- 
velopment of a ridge in ^(r) along a valley of | r/ (r ) | . 
The observation of it by NMR or /iSR experiments will 
form a direct evidence for the state II as well as for the 
presence of multi-order parameters. It is also quite in- 
teresting to perform the experiments in UPt3 where a 
lattice distortion has already been detected (26j. We fi- 
nally point out that the second-order transition between 
I<->II or I^->III is present for an arbitrary field direction 
in the basal plane. A glance on the functional (M may 
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FIG. 5. The vortex-lattice phase diagram for H || c as a 
function of \v\ and B/ H C 2 for ki = 1.2. The angle ip' L is j 
and for u>Q and v <0, respectively. 

lead to the conclusion that the transition I^-TII dis- 
appears for a low-symmetry direction, since the term 
|??i| 2 |?72| 2 yields those like vW^Wil 2 - However, it does 
persist as the transition (i) classified in the preceding 
paragraph. The hexagonal lattice has been checked to be 
stable in the high-field region, and the phase diagram for 
a small \v\ is qualitatively similar to Fig. 1. 

We finally present the results for H || c. Figure 5 shows 
the vortex lattice structure as a function of \v\ and B 
for K\ = 1.2. The square lattice is stabilized near H C 2 for 
small values of \u\, confirming Agterberg's result through 
a perturbation expansion with respect to v (k,\ = 1.2 cor- 
responds Agterberg's k^0.66 for v — Q) p2| . As the field 
is decreased, however, the lattice deforms into a rectan- 
gular one for \v\ < 0.17, followed by a further transition 
into the square and/or a distorted (i.e. 1; i9 ^ |, ^) 
lattice for \v\ < 0.1. The same calculation for k,\ = 2.6 re- 
veals that all the phase boundaries move rightward, with 
the distorted, square, and rectangular regions extending 
over 0<B/H c2 <l for |z^| < 0.004, 0.02 < \v\ < 0.09, and 
0.23 < \v\, respectively. With «i and \v\ small, the free en- 
ergies of these lattices are not much different from one an- 
other, as suggested by Agterberg's v-n diagram near H C 2 
[ p2| , and the present calculation reveals that there may 
also be field-dependent transformations among them. Al- 
though Riseman et al. |l7j| have reported an observation 
of the square lattice, there may exist field-dependent dis- 
tortion in the diffraction pattern. A detailed experiment 
on the field dependence may be worth carrying out. 
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